Modern NN Force Predictions for the Total ND Cross Section up 

to 300 MeV 



H. Witala*, H. Kamada*, A. Nogga , W. Glockle 
Institute for Theoretical Physics II, Ruhr- University Bochum, D-44 780 Bochum, Germany. 



Ch. Elster 

Institute of Nuclear and Particle Physics, and Department of Physics, 
, Ohio University, Athens, OH 45701 

G\ 

D. Hiiber 

Los Alamos National Laboratory, M.S.B283, Los Alamos, NM 87544 

(February 9, 2008) 



X 



Abstract 



> 

O _ 

For several modern nucleon-nucleon potentials state-of-the-art Faddeev calcu- 
. lations are carried out for the nd total cross section between 10 and 300 MeV 

projectile energy and compared to new high precision measurements. The 
agreement between theory and data is rather good, with exception at the 
higher energies where a 10% discrepancy builds up. In addition the conver- 
gence of the multiple scattering series incorporated in the Faddeev scheme 
is studied numerically with the result, that rescattering corrections remain 
important. Based on this multiple scattering series the high energy limit of 
the total nd cross section is also investigated analytically. In contrast to the 



5_i ' naive expectation that the total nd cross section is the sum of the np and nn 

total cross sections we find additional effects resulting from the rescattering 
processes, which have different signs and different behavior as function of the 
energy. A shadowing effect in the high energy limit only occurs for energies 
higher than 300 MeV. The expressions in the high energy limit have quali- 
tatively a similar behavior as the exactly calculated expressions, but can be 
expected to be valid quantitatively only at much higher energies. 
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I. INTRODUCTION 



For three-nucleon scattering the total neutron-deuteron (nd) cross section is the simplest 
observable, since it is integrated over the angular distribution in elastic nd scattering and all 
the angular and continuous energy distributions of the three-nucleon (3N) breakup process. 
If the theory does not agree with experiment, one has to expect that for some individual 
observables even stronger discrepancies might exist. On the other hand if there is agreement, 
possible discrepancies in some individual unpolarized differential cross sections have at least 
to average out. It is the aim of this article to compare precise new data between 10 and 300 
MeV neutron laboratory energy for the total nd cross section []J with fully converged Faddeev 
calculations based on the most modern nucleon-nucleon (NN) forces. The calculations are 
based on a strictly nonrelativistic treatment. Despite of this apparent restriction, we think 
that this comparison of theory and experiment will be an important benchmark result. 

A first presentation of our results appeared in [JIJ. Here we want to give a more detailed 
description and explicitly show that our theoretical results are stable under the exchange of 
one of the most modern NN potentials by another one. Detailed studies on the size of the 
contributions of the first few terms in the multiple scattering expansion to the total nd cross 
section are given. 

In order to obtain more insight into the behavior of the multiple scattering series for 
the total nd cross section we provide analytical high energy expansions for the first few 
terms. Compared to the most naive picture, which in the high energy limit would equate 
the total nd cross section with the sum of the total cross sections for neutron-proton (np) and 
neutron-neutron (nn) scattering, the rigorously calculated result is smaller for energies up 
to 300 MeV. Obviously shadowing effects can be expected. However, rescattering processes 
may a priori also enhance the total nd cross section. It is shown that rescattering processes 
of second order in the NN t-matrix lead to shadowing and antishadowing effects. In principle 
these features are known for the 3N scattering amplitude from studies in the framework of 
Glauber theory In contrast to this formulation we start from a multiple scattering 

series, which implies that it is an ordering according to powers of the NN t-matrix. The 
analytical steps leading to the high energy limit of the first leading terms are carried out 
in well defined integrals. There are no a priori assumptions involved about the scattering 
process, such as, e.g. diffraction type approximations. We observe that in the high energy 
limit a shadowing effect only occurs for energies larger than 300 MeV, below this energy the 
contributions second order in the NN t-matrix in the high energy limit enhance the sum of 
np and nn total cross sections. Furthermore, it is interesting to see that even at 300 MeV 
the high energy limit is quantitatively not yet reached. Nevertheless, the analytical studies 
make the underlying physics more transparent. An analogous investigation for the much 
simpler case of 3 bosons was performed in Ref. ||. 

The paper is organized as follows. In Sec. II we briefly describe the Faddeev framework, 
its multiple scattering expansion, the leading order terms in the NN t-matrix for obtaining 
the nd total cross section, and the high energy limit of the corresponding expressions. The 
derivation of the high energy limit is explicitly given in Appendices A and B. In Sec. Ill we 
present the results of our exact Faddeev calculations in comparison with the data and with 
the calculations in the high energy limit. In Sec. IV we study the effect of three-nucleon 
forces (3NF) on the total nd cross section. In addition, a simple estimate of relativistic 
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kinematic effects is given. We conclude in Sec. V. 



II. LEADING MULTIPLE SCATTERING TERMS FOR THE ND TOTAL CROSS 

SECTION AND HIGH ENERGY LIMIT 

We solve the Faddeev equations for 3N scattering in the form 

T|$) = fcP|$) + tPG T\$), (2.1) 

where t is the NN t-matrix, the operator P the sum of a cyclical and an anticyclical per- 
mutation of 3 objects and |<E>) the initial channel state composed of a deuteron and the 
momentum eigenstate of the projectile nucleon. The full breakup operator is 

U =(1 + P)T, (2.2) 

defining the physical meaning of the T-operator. The iteration of Eq. fl2.1|) together with 
Eq. ( |2.2p generates the multiple scattering series for the breakup process. 
Furthermore, the operator for elastic nd scattering is given by 

U = PGq 1 + PT. (2.3) 

Our aim is the evaluation and the understanding of the total nd cross section, which is given 
via the optical theorem 



=-(2vr) 3 -^Im($|f/|$>. (2.4) 
3<?o 



Here m is the nucleon mass and go is the asymptotic momentum of the projectile nucleon 
relative to the deuteron. Thus solving Eq. ( |2.1|) and using Eq. (|2.3|) yields the desired result 
in the form Eq. ( |2.4| ). 

In order to achieve insight how the total nd cross section is formed at high energies we 
consider the multiple scattering expansion of the operator U as given in Eq. ( |2.3j ). The first 
few terms constitute a power series expansion in the NN t-matrix 

U = PG 1 + PtP + PtG PtP + .... (2.5) 

At higher energies only a few terms in this expansion may be sufficient to describe the total 
nd cross section. We evaluate the different terms in the series numerically and study their 
contributions to the total cross section in Sec. III. 

In order to obtain analytical insight into the contributions to the total cross section we 
consider the imaginary part of the operator U as given in Eq. ( |2.5|) between the channel 
states |$) and obtain 

2iIm($|E/|$) = ($\U\$) - ($\U\$y 

= (<t>\P{t - t j )P\<S>) + (<5>\PtG PtP - Pt ] GlPt ] P\<S>) + ... (2.6) 

Since the first term in Eq. ( |2.5|) is real, it does not contribute to the total cross section. The 
second order term in the NN t-matrix can be rewritten as 
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($\PtG PtP - Pt ] GlPt ] P\§) = ($\P{t - rf)PG tP\$) 

- ($|P(t - tf)PG tP\$)* + (<S>\Pt^P{G - G* )tP\$), (2.7) 

which will be the starting point for extracting the high energy limit analytically. 

Let us first concentrate on the first order term of Eq. Q2.61 ) and work out the permutations. 
For the channel state |$) we choose 

|$) = = l^lqoh (2.8) 

and the NN t-matrix as given in the subsystem (23) 

t EE t X EE t 23 . (2.9) 

The index 1 is the convenient and usual notation |7|] to single out the pair (23). Similarily 
we shall denote t 2 = t 3 i and t 3 = tx2- The permutation operator P is given as 

P = + Pi A- (2.10) 

Next we define 

|$} 2 EE Px2P23\^)x 

m 3 = Pl3P2 3 mx (2.11) 

Applying P to the left and right one obtains for the first term on the right hand side of 
Eq. (p|) 

($|P(*-*t)P|$)= 3 ($|t-t t |$) 2 + 2 ($|t-* t |$> 2 

+ 3 ($|t_tt|$) 3+ 2 ($|t-tt|$) 3 . (2.12) 

An immediate consequence of Eq. (|2.11 ) and the antisymmetry of the deuteron state 
\<Pd) is 



|$> 3 = -P 23 |$) 2 . (2.13) 



Consequently one has 



3 ($|t-t t |$) 3 = 2 ($|P 23 (t-t t )P 2 3|$) 2 = 2 ($|t-t f |$) 2 

3 ($|t - ft|$) 2 = _ 2 ($|P 23 (t - tt)|$) 2 = 2 ($| t _ £ t|$) 3 . (2.14) 

Here we used the symmetry of t under exchange of particles 2 and 3. 

As final result we obtain for the first term on the right hand side of Eq. fl2.6|) , the first 
order term of the multiple scattering expansion for the nd total cross section 

($|P(t - t f )P|$) = 2 { 2 ($\t - t t |$) 2 + 2 {$\t - ^1^)3} 

= 2 2 ($|(t-tt)(i_P 23 )|$) 2 (2.15) 

It is a straightforward but tedious algebra to evaluate Eq. (|2.15 ) in the high energy 
limit. The derivation is sketched in Appendix A. Thus, the first order term in the multiple 
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scattering expansion gives in the high energy limit the following contribution to the total 
nd cross section 

4 1 iot = < t + C t - (2-16) 

The high energy limit is defined as the projectile momentum q being much larger than 
the typical momenta inside the deuteron. Thus Eq. ( |2.16|) shows that in the high energy 
limit the total cross section for nd scattering in first order in the NN t-matrix is indeed given 
as the sum of the total cross sections for np and nn scattering. This corresponds to the naive 
expectation, where the projectile nucleon scatters independently from both constituents of 
the deuteron. 

Next we study the rescattering processes of second order in t as given in Eq. ( f2.7| ). 
Working out the permutations the first term of the right hand side of Eq. (|2.7Q can be 
written as 

($|P(t - tf)PG tP\$) = 2 { 2 ($|(t - t^Gohd®), + |$> 2 ) 

+ 2{ 2 (m-J)G t 2 m i + \$) 3 ). (2.17) 



Similarly the remaining term in Eq. (|2.7|) yields 

($\PtfP(G - G* )tP\$) = -4vu { 2 ($\tf6(E - F Q )f 3 (|$) 1 + |$) 2 ) 

+ 2 (^8(E-H )t 9 (\9) 1 + \9) a )}. (2.18) 

The analytic evaluation of the high energy limit is even more tedious and is sketched in 
Appendix B. The leading contribution of the terms given in Eq. ( |2.7| ), second order in t in 
the multiple scattering expansion to the total nd cross section, is given by 



1 / f°° 
■7-Y?~ l l dr(pi(r)(pi. 

[2 c(iiz, i, — i) c(m',i,-i) 

, 1.4m. 9 . Nfi/ 1 /T ,3 11 ,3 11 . 

--(Re(--qo--n P \t\- qo --np) a ) 

1 /T , 3 11 . .3 1 1 , N2 

+ Im \~i q ° ~ 22 nP '*'l qo 2 ~~ 2 nV ' a > 
!/t, / 3 11 . ,3 1 1 . , 2 
Re ^ 4 q ° ~ 2 2 nP ' 'Z q °2 ~ 2 nP ' a > 

„ ,3 11 . ,3 11 . „ ,3 11 . .3 11 , 
+Re(-qo--np|t|-qo--np) a Re(-q --nn|t|-qo--nn) a 

„ ,3 1 1 ,.3 1 1 , „ ,3 1 1 . .3 1 1 . N , 
+^(4^0- - -nn\t\-q - - -nn) a Re(-q - - -np|*|-q - - -np)J} 

+C(ll/,00) c(iir,oo) 

, l,4m. 9/ ...1 /T . 3 1 1 .3 1 1 . .9 
{ -2W (2X) ( 2 (Im <-4 q °2 - 2" p|t| 4 q °2 - 2" P> J 

1/t, / 3 1 1 , , 3 1 1 ^2 

~2~( Re \ 4 q °2 ~ 2 nP ' '4 q °2 ~ 2 nP ' a > 
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T , 3 11 . ,3 11 . T . 3 11 , , 3 1 1 
+lm{--q --np\t\-q --np) a lm{--q - ^np\t\-q - - -np) a 

T /3 1 1 ,3 11 . T ,3 1 1 , ,3 11 . 
- Im (^qo2 - ^W^o _ 22^ a Im M q °2 ~ 2 nn ^4 q ° ~ 22 nn ' a 

„ .3 11 . .3 11 , „ .3 1 1 . .3 1 1 " 
+Re(-q 2 2 n Pl t l4qo2 2 n Wa Re (4qo2 - -^nn^-^- ~ 2 Wl ' a 

„ ,3 11 . .3 11 . „ .3 1 1 ,3 1 1 . 
+^\^o^nn\t\-qo--nn) a Re{-q Q - - -np\t\-q - - -np) a 

n , 3 11 M 3 11 .„ , 3 11 , , 3 1 1 
-Re(--q --np|t|-q --np) Re(--q - -~np\t\~q - - -np) a 

n ,3 1 1 . ,3 " 1 1 . „ ,3 1 1 . ,3 11 . N11 
+ Re (^qo2 - ^^q^ _ 22 n ^ a Re U q °2 ~ 2 nn ' 'l q ° ~~ 22 nn ' a) ^ 



1 / f°° 

drtpi{r)<p v \ 



r) x 



4vry Jo 

[2 c(m,i,-i) c(iu',i,-i) 

I — (rr tot ( — —)rr tot ( — — ) 4- rr tot (— — — )rr tot (— — — 

\ 2 v np V2 2 ' nn 22 n 2 2 2 2 

+c(n/,oo) c(nr,oo) 



r 1 / tot/ 1 1 \ tot / 1 1\ 1 1 \ ^tot / 1 

l^npl^ nn ^2 ~~ r + «p^2 ~~ 2^ nn ^22^ 

txf + afi (2.19) 



Here the sum over 1, 1' denote the S and D-wave contributions of the deuteron wave functions 
(fi(r). Further occur the imaginary and real parts of the forward and backward two-nucleon 
scattering amplitudes (nn or np) with specified spin magnetic quantum numbers in the 
initial and final states. The index a denotes antisymmetrization without the factor The 
terms of second order in t give rise to Eq. ( |2.19| ), a quite complicated expression exhibiting 
the interferences due to spin and isospin degrees of freedom. The part a) related to the 
first square bracket in Eq. Q2.19p contains real and imaginary parts of the NN scattering 
amplitudes and has positive as well as negative contributions. Only a numerical evaluation 
of these terms as function of the scattering energy can give insight about the size and energy 

(2) 

dependence of these terms. The part cr}/ related to the second square bracket contains due 

(2) 

the optical theorem a product of two NN total cross sections. That part a}/ is negative and 
reduces the term of first order in the multiple scattering expansions, the sum of the two NN 
total cross sections. This is naturally called shadowing effect. 

Neglecting spin and isospin dependencies and consequently the D-wave admixture in the 
deuteron, and further neglecting the backward amplitudes that expression given in Eq. ( |2.19|) 
reduces to the much simpler ones which are given in || for 3 bosons. 

III. DISCUSSION OF THE RESULTS OF THE FADDEEV CALCULATIONS 
AND OF THE HIGH ENERGY LIMIT 

For calculating the nd total cross section we employ the most recent NN potentials, 
namely CD-Bonn [[J, AV18 [1(| ,Nijm I and II [[IT]]. Those potentials are optimally fitted to 
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the Nijmegen data basis up to 350 MeV nucleoli laboratory energy. The Faddeev equation 
Eq. (2.1) is solved in momentum space in a partial wave decomposed form. For a description 
of the technical details of the calculations we refer to Ref. H. In Fig. |I] we compare the 
calculated nd total cross section based on the CD-Bonn potential with the data. The 
experimental data are obtained by adding the separately measured values of the hydrogen 
cross section and the deuterium-hydrogen cross section difference 0. The figure shows 
that the theoretical calculation describes the data very well at the lower energies but falls 
below the data at higher energies. The calculation begins to underestimate the data around 
100 MeV by about 4%, and this discrepancy increases to about 11% at 300 MeV. 

The calculation presented in Fig. [1| employs in the two-nucleon subsystem angular mo- 
menta up to jmax=4:. This is sufficient if we are satisfied with a calculational accuracy of 
1%, which corresponds to the size of the experimental error. In order to demonstrate the 
dependence of the numerical accuracy of our calculations, we show in Table | for a few 
energies in the relevant energy regime the convergence of the total cross section as func- 
tion of jmaxi the maximum angular momentum of the two-nucleon subsystem. In the 3N 
calculations two-nucleon angular momenta larger than j max are put to zero. We see that 
jmax = 4 is sufficient, if we are satisfied with an accuracy of 1%. Next we demonstrate the 
stability of our theoretical result under the exchange of the NN potentials. This is shown 
in Table [II]. Clearly within an accuracy of 1% the predictions of those four essentially phase 
equivalent potentials agree with each other. From this we can conclude that the deviation 
of the calculation from the data at higher energies, which is shown in Fig. [I] is independent 
from the NN interaction employed. 

In order to gain some insight how the total nd cross section is build up, the individual 
contributions of first, second and third order in the multiple scattering expansion of the 
elastic amplitude are compared. In Fig. |2| we present the total cross section calculated in 
first order, successively add the contributions of the second and third order and compare 
those results with the full calculation. The figure shows that the first order in the multiple 
scattering expansion is larger than the full result at all energies. The second order contribu- 
tion enhances the cross section below about 130 MeV and decreases it at the higher energies. 
We also see that at energies below about 200 MeV the third order term is significantly larger 
than the second order term, and thus that rescattering of higher order is very important. 
Only above 200 MeV the third order contribution becomes negligible. Here the small cor- 
rection of the second order rescattering process becomes sufficient to describe together with 
the first order contribution the full solution of the Faddeev equation. Around 300 MeV the 
first two terms in the multiple scattering expansion of the scattering amplitude are sufficient 
to describe the total cross section. 

As a next step we consider the high energy limit of the first and second order terms 
of the multiple scattering expansion as derived in Sec. II. The first order term results in 
a sum of the total cross section for np and nn scattering. The high energy limit of the 
second order term is given in Eq. Q2.19 ). In this equation we already indicate that this term 
consists of two parts, one which consists of products of real and imaginary parts of the NN 
t-matrix, and one which is due to the optical theorem proportional to products of NN total 
cross sections. The latter term is negative and can be considered as shadowing effect. The 
contributions of the different terms of the high energy limit are shown in Fig. We show 
the individual contributions up to a laboratory energy of 800 MeV in order to get a better 
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insight into their energy dependence. The contribution of the sum of a np + a nn is by far 

(2) 

dominant. The calculation shows that a) is always positive and thus enhances the total 

(2) 

cross section. Thus this term could be viewed as anti-shadowing effect. The term <j\{ is 

(2) (2) 

always negative. We see that a\ falls off faster as function of the energy than a)/. This 
interplay of the different contributions to the high energy limit of the second order term 
in the multiple scattering expansion leads to an enhancement of the total nd cross section 
below about 400 MeV and a weakening above compared to the sum of the np and nn total 
cross sections. Thus only above 400 MeV the second order term in the elastic amplitude 
leads in the high energy limit to a shadowing effect as it is familiar from older estimates 
within the framework of Glauber theory. This is shown in Fig. FJ. The dashed line gives the 
contribution of the sum of a np + a nn . To this sum the term a\ is added and given by the 

(2) 

dotted line and enhances the total cross section. Then we add in addition a n leading to 
the solid line. 

It is now interesting to compare those results obtained from the first two terms of the 
multiple scattering expansion in the high energy limit with the corresponding contributions 
calculated exactly with the Faddeev framework. This comparison is summarized in Table [TTI] . 
The second column gives the exact Faddeev result as reference. It should be noted that all 
results given in Table [Til] are calculated using j max =4 in the two-nucleon subsystem. Starting 
at 100 MeV we compare the exact first order result to its high energy limit, the sum of the 
np and nn total cross sections, as given in Eq. (|2.16| ). The table shows that this sum is larger 
than the exact first order term for 140 MeV and higher. In fact, the difference increases 
with increasing energy. In the case of three bosons interacting via Malfliet-Tjon type forces 
(acting in all partial waves) this difference vanishes around 300 MeV and the asymptotic 
expression acquires its validity here ||. Apparently, the three nucleon system with realistic 
forces behaves differently. It would be very interesting to find out at which energy the first 
order term will reach the asymptotic form. This interest is however restricted to the context 
of a mathematical model of potential scattering, since the standard NN forces loose their 
physical meaning above pion-production threshold. 

Next we compare the exact second order result to the expression of the high energy 
limit as given in Eq. ( |2.19|) . First, the exact second order result of the Faddeev formalism 
is a small correction to the first order result. At lower energies it is additive and changes 
sign around 150 MeV, which then leads to a reduction of the cross section. As can be 
seen in Fig. fj, the second order term in the high energy limit exhibits this behavior only 
around 400 MeV. Comparing the sum of first and second order terms of the exact Faddeev 
calculation to the corresponding term in the high energy limit (column 9 which is the sum 
of the terms in columns 6-8 in Table [IT]] ) shows that the high energy limit is always larger in 
the energy region under consideration. At 300 MeV the difference for the combined results 
in those orders is still 8%, and thus the high energy limit is not yet quantitatively valid. 

As remark we would like to add that the asymptotic form requires higher NN force 
components than j max =4. For the study presented in Table [nj we kept the angular momenta 
in the two-nucleon system fixed at this value in order to compare with the Faddeev results. 
The two-nucleon t matrices underlying the calculations presented in Fig. however, employ 
up to j m ax = 15 in order to obtain a converged result at 800 MeV. In the case of the high 
energy limit this can be easily achieved, since only the NN on-shell amplitudes enter. 
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IV. ESTIMATE OF THREE-NUCLEON FORCE AND RELATIVISTIC 

CORRECTIONS 



In this section we want to study the effect of three nucleon forces (3NF) on the total nd 
cross section and give a very simple estimate on the size of relativistic effects. The inclusion 
of 3NF in 3N continuum Faddeev calculations has been formulated in Ref. [|1^], superseding 
a previous formulation given in Ref. || . The newer formulation is formally more elegant and 
numerically more efficient and is used in our calculations. The inclusion of 3NF generalizes 
Eq. to 

T|$) = iP|$) + (1 + tG )vi l \l + P)|$) 

+ tPG T\$) + (1 + tG )V 4 (1) (l + P)G T\$), (4.1) 

where V4 is that part of the 3NF, which is symmetric under exchange of particles 2 and 3. 
The total 3NF then has the form 

V 4 = vi 1] + V} 2) + V} 3 \ (4.2) 



and is totally symmetric. The full breakup operator as given in Eq. fl2.2| ) keeps its form, 
whereas the operator for elastic nd scattering changes to 

U = PG 1 + V} 1} {1 + P) + PT + V} 1] {1 + P)G T. (4.3) 



The optical theorem given in Eq. (|2.4j) stays of course valid. 

The numerical evaluation of Eq. ( |4.1|) is quite demanding, since in contrast to the 3N 



bound state one needs the partial wave projected momentum space representation of V4 
now for both parities and, in addition, for total 3N angular momenta larger than J = 1/2. 
At a projectile energy E™ ah = 200 MeV we checked that total angular momenta up to 
J = 9/2 were needed. The calculation of the three-nucleon force matrix elements for such 
high angular momenta became possible due to the new partial wave decomposition for the 
three-nucleon force introduced in [T3|. Of course, NN forces contribute significantly up to 



J = 25/2 to 3N scattering states. In our calculations we employ the Tucson- Melbourne 
force [HJ], which has has been adjusted individually to the 3 H binding energy for each of 
the different modern NN forces |T3|. The results presented in this paper are based on the 



CD-Bonn potential and the corresponding parameters of the Tucson-Melbourne force. 
In a previous study we found that 3N force effects are responsible for filling the minima in 



the angular distribution of elastic nd scattering, especially at higher energies [1J|. Though 
the effect is seen only in the minima of the differential cross sections, and thus is small 
in magnitude, one might expect traces thereof in the integrated quantity, the total cross 
section. Corresponding effect to those in the elastic angular distribution might also occur in 
3N breakup processes, over which is also integrated when calculating the total cross section. 
This is indeed the shown in Fig. [5| The effect of the 3N force enhances the total 

nd cross section calculated with NN forces only by about 4%. It is interesting to note that 
this enhancement is almost independent of the energy when considering the energy regime 
between 100 and 300 MeV. It is also obvious from the figure that the 3N force effects included 
here are not strong enough to bridge the gap between the calculation based on 2N forces only 
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and the data. However, one should keep in mind that the Tucson-Melbourne 3NF model is 
just one model, and theory as well as application of 3N forces is just at the beginning. 

A second effect, which can be expected to appear especially at higher energies is an effect 
due to relativity. A generally accepted framework for carrying out relativistic 3N scattering 
calculations does not yet exist. Thus we restrict ourselves to a very simple estimate of a 
relativistic kinematic effect. Due to the optical theorem the total cross section is given as 

o£ X ±r1m{*\U\Q), (4.4) 

where \j\ is the incoming current density. If we neglect possible changes of Im($|[/|$) in 
going from a nonrelativistic to a relativistic formulation, the ratio of the relativistic to the 
nonrelativistic total cross section is simply given by 

u tot,rel \j\nr E n Ed , m n m d , . 

cCt,nr ~ |7U ~ p r n el (En + E d ) / p^(m n + m d y 

Here, E n and E& are the relativistic kinetic energies of the neutron and the deuteron in the 
cm. system, and p r ^ 1 is the relativistic cm. momentum. The other terms are obviously 
the nonrelativistic approximations thereof. This ratio given in Eq. (4.5) is easily evaluated 



and leads to an increase of the total nd cross section, namely an increase of about 3% at 
100 MeV and about 7% at 250 MeV. Calculations based on 2N forces only, which include 
this relativistic effect, are shown at various energies in Fig. [|. Again, these additional 
effects enhance the total cross section and bring the calculation closer to the data. It is also 
clear that this effect does not constitute a relativistic theory, since the forward scattering 
amplitude is still calculated entirely nonrelativistically. However, both effects, the one caused 
by inclusion of a 3NF model and the simple estimate of a relativistic effect, are roughly 
similar in magnitude, and when added up would come close to the data in the higher energy 
regime. It should be also noted that at lower energies, the added effects would overpredict 
the data. 



V. SUMMARY 

In view of new precise experimental information on the total nd cross section fully con- 
verged Faddeev calculations based on the most modern NN forces are carried out at projec- 
tile energies between 10 and 300 MeV. For the Faddeev calculations a strictly nonrelativistic 
treatment is employed. The calculations show that the results for the total nd cross section 
do not depend on the choice among the most recent phase-equivalent potentials. 

In order to obtain more insight into the behavior of the multiple scattering series for 
the nd total cross section, we study its convergence within the Faddeev framework. Below 
about 200 MeV projectile energy rescattering of higher order is very important, as can be 
concluded from the fact that here the third order contribution of the multiple scattering 
series is significantly larger in magnitude than the second order. Only around 300 MeV the 
first two terms in the multiple scattering expansion are sufficient to describe the total nd 
cross section. 

In order to obtain analytical insight we investigate the first terms of the multiple scat- 
tering series for the forward elastic nd scattering amplitude resulting from the Faddeev 
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equations in the high energy limit. Although the treatment is purely nonrelativistic and 
enters far into the region where relativity is expected to be important, we think it is inter- 
esting to know the asymptotic behavior at high energies for pure potential models without 
absorption. Absorption processes (particle production) occurring in a relativistic context 
will change the results presented here. In accordance with the naive expectation, we extract 
from the first order term in the multiple scattering series the sum of the total cross sections 
for np and nn scattering, and from the second order term a shadowing effect proportional to 
products of the total NN cross sections. The shadowing contribution is negative and reduces 
the total nd cross section. However, we also find positive terms proportional to products of 
real and imaginary parts of spin dependent NN forward and backward scattering amplitudes. 
Those terms provide an overall enhancement of the nd total cross section, but vanish faster 
with energy than the shadowing term. At projectile energies above 400 MeV the shadowing 
is the dominant second order effect. 

A comparison of the theoretical calculations based on a nonrelativistic Faddeev frame- 
work with two-nucleon forces only and the experimental observables exhibits a discrepancy 
with respect to the nd data, which starts around 100 MeV with a few percent and reaches 
about 10% at 300 MeV. Possible corrections can be due to either three-nucleon forces or 
relativistic effects or both . We calculated the effect of the Tucson-Melbourne 3N force on 
the total nd cross section and find that in the energy regime between 100 and 300 MeV 
this 3NF model provides an overall enhancement of the total cross section of ~4%. This 
is consistent with previous studies, which found that the minima of the elastic scattering 
angular distribution are filled in when this 3NF is taken into consideration. In the total 
cross section this effect is not large enough to bring the present calculations close to the 
data, but the trend is in the right direction. 

A second effect, which can be expected to appear especially at higher energies is due to 
relativity. We carry out a simple estimate of effects due to relativistic kinematics , and find 
that these corrections enhance the total nd cross section about 7% at 250 MeV, while being 
much smaller at lower energies. Again, those corrections go into the right direction with 
respect to the data. However, we do not suggest that our estimate is the complete solution 
of the problem, since in our calculation the forward scattering amplitude is still calculated 
entirely nonrelativistically. 

Both our estimates of corrections to our nonrelativistic Faddeev calculations based on 2N 
forces suggest, that effects due to 3N forces as well as relativistic effects become non negligible 
at higher energies. This calls for a strong effort to progress theoretical developments of a 
theory of 3N forces as well as a relativistic framework for 3N scattering. 
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APPENDIX A: THE HIGH ENERGY LIMIT OF THE FIRST ORDER TERM IN 

THE NN T-MATRIX 



To evaluate the first order term of the multiple scattering expansion for the total cross 
section in the high energy limit we consider 

Xl EE 2 ($ | (*-*t) (1-P 23 ) |$> 2 . (Al) 

The completeness relation in the three-nucleon space is given by 
1= E E [d 3 P d 3 q 

| p m 2 m 3 v 2 v 3 > |q m\V\ > < p m 2 m 3 z/ 2 z/ 3 | < q m\V\ | (A2) 

Here p and q are standard Jacobi momenta, the spin magnetic quantum numbers 
and fi the corresponding isospin quantum numbers of the three nucleons. In that notation 
the channel state |<3>) 2 reads 

1$ >2 = \<Pdm d Vd > |qo m N u N > . (A3) 
We also need a change of basis 

2 < p m^mi v 3 V\ q m 2 z/ 2 |p m' 2 m' 3 z/ 2 z/ 3 q ' m'^ > 1 

1 1 3 3 

= <Kp ' - - q' - q) 5(p + q' + - q) JJ <5 mjm > J] <M ( A4 ) 

2 2 i=i i=i 

and the explicit expression for the deuteron state 
< p m 3 mi u 3 ui\ip d > = 

7k— $vi-v3 E (l^,mim d -mi) 

Y lmi ip) 2 1 '" 23 mi md ~~ m ^ ( ^ A5 ' ) 

As an immediate step one finds after a suitable substitution of integration variables 
2 < $|i - t f |$ > 2 = E EEE / ^<V<i|-z m 3 miz/ 3 z/i > 

3 1 tl 3 1 

< ^qo + -zm N m 3 u N u 3 \t - t 1 1 -q + -zm N m 3 v N v 3 > 

< —z.m! 3 miv 3 vi\(p d > 

If the projectile momentum | qo is much larger than typical momenta z occurring in the 
deuteron state, one can take the NN t-matrices out of the integral and obtains 



15 



2 < <E>|£ - t f |$ >2 



I go— >oo 



3 3 
H Him < ^qo^iY^3^^3^ - t ] \-^m N m' 3 v N v' 3 > x 



m i m 3 v i v 3 ^ 



y <i 3 2; < (pd\ — 7,m 3 m\v 3 v\ > < —zm 3 miv 3 isi\(pd > . (A6) 



Here the energy argument of the t matrices is given by 



3 o ^ , 1 o 9\ 

= £rf+ 4^ q °-4^ ( 4 q °- qo - Z + Z } 

1,3, o 3 3p 
= e d + — (-q ) +— q -z-— z 

m 4 4m 4m 

>tr I ( i_ qo)2 . (A7) 

m 4m 

Thus in the high energy limit we only encounter on-shell NN t-matrices. 
Using Eq . (|A5|) we find 



3 3 
2 <$|t-t t |$>2 — > X] < j^omNmsUN,-^ -t *\-q.om N m 3 v N ,-vi > 

rnims v\ 
\ poo 

o / d P P 2 vKp) C(lll, m d - m 3 - mi, m 3 + mi) 2 
2 T io 

C(~l,m 8 m 1 ) a . (A8) 

The next step is to introduce properly antisymmetrized NN states in the NN t-matrix. This 
comes automatically of course, since the form X\ contains the operation (1 — P 23 )- Taking 
this into account we find altogether for the expression Eq. (|2.15|) 

3 3 
< $\P(t - P )P|$ > — > na < -q m N m 3 u N , -v x \t - P\-q m N m 3 v N , -v x > na 



poo 

V / dp p 2 (pf(p)C(lll, m d - m 3 - mi, m 3 + mi) 5 
, Jo 



x 

i 

C(iil,m 3 mi) 2 . (A9) 
Here the antisymmetrized and normalized free NN state is given by 

|pm 2 m 3 z/2 v 3 > na = -j=(l - P 23 )\pm 2 m 3 v 2 v 3 > . (A10) 

Finally we use the relation 

3 m x — , ^ — , r — .3 . , I,, tj\ t\\ 3 



t - tt 



27li 4 g °Y J2 J dpV\-qQpm' 2 rn' 3 v' 2 v 3 > (+) (+) < -q pm' 2 m' 3 u 2 u' 3 \V, (All) 



rn:,iii,, v\v. 
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where the NN potential V is applied to the two-nucleon scattering states | . . Further- 
more, the definition of the total NN cross section initiated by fixed magnetic spin and isospin 
quantum numbers is given by 



m 2 m 3 V 2 U 3 

(A12) 

This inserted into Eq. ( |A9| ) and applying the relation given in Eq. (|2.4j ) to the total nd cross 
section we arrive at 

ffSK,vs)= 12 12 a l°N,-vi( m N, m 3) 

1711)713 V\ 

roc 1 1 



X 



/-oo J I 

/ rf P V 2( fi (p)C(Zll, m d -m 3 - mi, m 3 + mi) 2 C(--, m 3 , mi) 2 . (A13) 
l Jo 2 2 



This is the contribution of the term Eq. Q2.15| ) to the total nci cross section for initially 
polarized particles and a nucleon species of type v N . Averaging over the initial state po- 
larizations one can perform the summations over m d and analytically and using the 
normalization conditions for the s- and d-wave parts of the deuteron wave function one ends 
up for a neutron induced process with 

tot „tot I „tot / A 1 A \ 

°nd = °np + CT nn- ( A14 ) 



APPENDIX B: THE HIGH ENERGY LIMITS OF THE SECOND ORDER 

TERMS IN THE NN T-MATRIX 

To evaluate the second order term of the multiple scattering expansion for the total nd 
cross section in the high energy limit we consider 

X 2 = 2 (m-t^) G t 3 |$)i (Bl) 
Corresponding steps to the ones carried out in Appendix A lead to 

X 2 = Yl 12 12 d3 z d 3z ' < Vd\ - zmgmiZAjZ/i 



> 



X 



i 2 u 2 ^3 m^V2 

< |Qo + | z m\m 3 v\v 3 \t - ^||q - ~z - z 'm^m^u^ > 
-\e d \ - ±(z 2 + z' 2 + z-z') + JL(z + z • q + xe 

x < ~^qo + z + 'mim' 2 vn/ 2 \ t |^q + ^z 'm^m^X' > 

< -z 'm^m^Xlv^ > • (B2) 

Again we assume that | q is much larger than the momenta z and z' in the deuteron state 
and find 



17 



X 2 — ► < j<lornlm 3 v%v 3 \t - t^-q rn! 2 m' 3 i>' 2 i>' 3 > 

x < — q m 1 m 2 z/ 1 z/ 2 | r |-q m 1 m 2 z/ 1 z/ 2 > -i, (,B3) 



4 

with 



J = / d A V I d 3 p' , -r. r.^^ ^yrr.r^ _ (fi4) 



< <fd\ 


— pm 3 m 1 i/3^ 1 >< — p 'm^m'gZ/^Vg 








-£(p a + p' 2 + p-p') + 5fe(p + p') 


• qo + «e 



For the extraction of the asymptotic limit as §|qo| — » oo we refer to a previous study 
carried out for 3 bosons. We use Eq. fl^5| ) and use the deuteron wave function given in 
configuration space. After some algebra one arrives at 



3g 

1 1 



x ^C(lll,0m d )(--l,m 3 m d -m 3 )V2l + T 
i 1 1 

x i l ~ J dr cpi(r)(pi<(r) 5 mi , md - m3 5 m ' 3>md - m > 2 . (B5) 

The second term in Eq. ( |2.17| ) can be evaluated analogously and we end up with the inter- 
mediate result 

($|P(t-tt)P Go tP|$) = -^(2vr) 2 £ £ £ 

mafg m »^' 

< j^. m 2 m 3 u 2 iy 3 \t - P\-q m 2 m d - m 2 , v 2 , -v 2 > na 

3 3 ii 

< --q m d - m 3 ,m 2 , -v 3 , v' 2 \t\-^m\m 2 v\v 2 > na (-) a -Ms (— )a-»^ 

5] V2/ + 1 C(Zll,0m d ) (~-l,m 3 m d -m 3 ) 

,11 



v / 27 7 + TC(/ / ll,0m d ) ( l, m > d -m 2 ') 



i'-' / dr W (r)^(r). (B6) 

JO 

We encounter only forward or backward NN scattering amplitudes, which induces connec- 
tions between the magnetic spin quantum numbers occurring in the t-matrices (put % = z). 
Further we sum over the initial spin magnetic quantum numbers m d and m^. As an example 
of these summations we show the result 

{(®\P(t - G tP\$) - complex conjugate} = 



m N m d 



roo 

V\/2TTTVv / 27 7 TTz'-W dr 

l V Jo 
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2C(ai,01)C(ril,01){< Ar (Hnp) +<7j&(iinn) + - inp)oj&(i - ^n)} 

+C(ni ) 00)C(ril,(X)){^(~np)a^ r (i - ~nn) +^(~ - ^^(—nri)} 



5m 
3g, 



i— (2vr) 2 V V2Z + 1 V; VWTl % l ~ l / <pi{r)(p v {r) 



2C(ui,oi)c(m,oi) 

3 11 . .3 11 N2 , t 3 1 1 . . 3 1 1 , 21 

{(Im < -q --np|t|-q --pn > na ) + (Jm < -qo^ - T^V I < I ^qo^ - ^ >na ) ' 

+c(m, oo)c(m, oo) 

3 1 1 . ,3 1 1 .. T 3 1 1 . .3 1 1 

{2(Im < -q --np\t\-q --pn > na )(Im < -q - - -pn\t\-q - - -np > na ) 

3 113 11 
+(Im < ~q - - -np\t\ - -q --pn > na f (B7) 

3 1 1 . . 3 11 wr 3 1 1 . . 3 11 

-2(Im < -q - - -np\t\ - -^--np > na ){Im < -q„- - -nn\t\ - -q --nn > na )\. 



The last term in Eq. Q2.7D can be treated analogously to Eq. (|Bl|) resulting in 
(^P(G -G* )tP\^)^-8m E E E 

m\mzv\V3, m 2 i m s r v 2 i m 2 iiv 2 ii 
^3 o I -tl / / / / . 

3 3 
x < -^qomim 2 ^i^ 2 |*|^qom?m 2 V°i/2 > na I, (B8) 



where 



I — J d 3 p J d 3 p ' < (p d \ — \>m 3 m\v 3 v\ >< — p m^m'^v^Wd > 

<K"M - -(p 2 + p ' 2 + p ' • p) + ^-q • (p + p '))• 
m 2m 

It is easily shown that in the high energy limit 

I= — I. (B9) 

ITT 

It remains to perform the summations over the spin- and isospin quantum numbers. Then 
using the connection to the total nd cross section, Eq. ( |2.4| ), the lengthy expression Eq. (|2.19 ) 
results. 
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TABLES 



^lab V ivlt; v ) 


<d (mb) 




jmax — 1 


jmax — 2 


jmax — 3 


jmax — 4 


jmax — 5 


10.0 


954.9 


1038.0 


1035.2 


1037.3 


1036.7 


60.0 


141.0 


177.8 


177.0 


177.4 


177.3 


140.0 


43.0 


71.6 


73.2 


74.7 


74.6 


200.0 


31.7 


55.8 


58.3 


60.5 


60.4 


260.0 


29.7 


49.4 


52.2 


54.7 


54.7 


300.0 


29.9 


46.9 


49.7 


52.3 


52.4 



TABLE I. The convergence of the calculated total nd cross section with increasing j rnax at 
selected energies. As NN interaction the CD-Bonn potential has been used. 



Eiab (MeV) 


< d (mb) 




AV18 


CD-Bonn 


Nijml 


Nijmll 


10.0 


1039.0 


1035.2 


1038.0 


1038.3 


100.0 


97.7 


99.3 


98.2 


97.5 


140.0 


72.2 


73.2 


72.5 


71.8 


200.0 


57.8 


58.3 


57.9 


57.4 


300.0 


49.5 


49.7 


49.6 


49.3 



TABLE II. Stability of the calculated total nd cross section a l °^ under the exchange of the 
NN potentials. The calculations were performed with jmax — 3. 



20 



^lab l ivie v ) 








tot 


mb) 














Faddeev Calculat 


ion 


High Energy Expansion 




full FC 


Lord. 


2.ord. 


l.+2.ord. 


&np + CT nn 


°? 


a ll 


sum 


100. 


100.95 


108.27 


1.91 


110.18 


105.60 


20.62 


-6.54 




119.67 


140. 


74.85 


78.62 


0.08 


78.69 


80.23 


12.20 


-3.96 




88.45 


200. 


60.55 


62.74 


-0.86 


61.89 


65.86 


6.28 


-2.80 




69.34 


300. 


52.30 


54.05 


-1.52 


52.53 


58.14 


2.38 


-2.28 




58.24 



TABLE III. Comparison of Faddeev results with the asymptotic expansion ones. The calcula- 
te) 

tions were performed throughout with the CD-Bonn potential and j ma x = 4. The terms a\ and 

(2) 

a \i are from Eq. (2.19) . For more explanation of the different terms see Sec. II and III. 
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FIGURES 



FIG. 1. Comparison of the Faddeev calculation for the total nd cross section based on the 
CD-Bonn || potential with data Q. The experimental values for the total cross section for deu- 
terium were obtained by adding the separately measured values of the hydrogen cross section and 
the deuterium- hydrogen cross section difference (for details see Ref. Error bars are omitted, 
since they are smaller than the dot size. 

FIG. 2. The contributions to the total nd cross section of the different orders in the multiple 
scattering expansion of the Faddeev amplitude. Successive orders are added to the first order term 
and then compared with the full Faddeev calculation. The potential employed is the CD-Bonn 
model. 

FIG. 3. The contributions of the different orders in the NN t-matrix in the high energy limit 
to the total nd cross section. The solid line shows the sum of the np and nn total cross sections. 
The positive contribution a\ (dashed line) drops faster in energy than the negative contribution 

(2) 

a j j (its magnitude is shown as dotted line) of the second order term of the multiple scattering. 
For details see Sec. III. 

FIG. 4. The contributions to the total nd cross section in the high energy limit. The dashed 
line shows the sum of the np and nn total cross sections. Successively added to this is the positive 

(2} (2) 

contribution a\ (dotted line) and the negative contribution crj/ (solid line). For details see 
Sec. III. 

FIG. 5. Corrections to the Faddeev calculation based on a nonrelativistic Hamiltonian using 
strictly two-nucleon forces (solid line) for the total nd cross section. The open squares show 
calculations at various energies, where the Tucson-Melbourne 3N force has been included. The 
open triangles show at various energies calculations based on two-nucleon forces corrected with a 
relativistic kinematic (for details see Sec. IV). The dots describe the data All calculations are 
based on the CD-Bonn potential. 
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